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Abstract 

We study dynamics of two coupled periodically driven oscillators. The 
internal motion is separated off exactly to yield a nonlinear fourth-order 
equation describing inner dynamics. Periodic steady-state solutions of 
the fourth-order equation are determined within the Krylov-Bogoliubov- 
Mitropolsky approach - we compute the amplitude profiles, which from 
mathematical point of view are algebraic curves. 

In the present paper we investigate metamorphoses of amplitude pro- 
files induced by changes of control parameters near singular points of these 
curves. It follows that dynamics changes qualitatively in the neighbour- 
hood of a singular point. 

1 Introduction 

In this work we study dynamics of two coupled oscillators, one of which is driven 
by an external periodic force. Equations governing dynamics of such system are 
of form: 



TJl\X\ - V\ (±i) - Ri (xi) + V 2 {%2 ~ ±l) + R 2 {%2 - X\) = / cos (ujt) 
m 2 x 2 - V 2 (x 2 - X\) - R 2 (x 2 — X\) = 



(1) 



where R\ , Vi and R 2 , V 2 are nonlinear elastic restoring force and nonlinear force 
of internal friction for mass mi and mass m 2 , respectively. Dynamic vibration 
absorber, consisting of a (generally small) mass m 2 , attached to the primary 
vibrating system of (typically larger) mass mi is a generic mechanical model 
described by © [HCI|. 
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We shall consider a special case: 

Ri (xi) = -aixi, Vi (x x ) = -v x x\. (2) 

Dynamics of coupled periodically driven oscillators is very complicated [3J 
HI [5J [S] . Starting from equations (fT|), ^ we derived the exact fourth-order 
nonlinear equation for internal motion as well as approximate second-order effec- 
tive equation [HI OH] (this approximation performs well for ^ <C 1). Applying 
the Krylov-Bogoliubov-Mitropolsky (KBM) method to the effective equation 
we have computed and studied the corresponding nonlinear resonances. More 
exactly, we investigated the amplitude profiles (resonance curves) A (w), i.e. de- 
pendence of the amplitude on the frequency u, given implicitly by the KBM 
method. Metamorphoses of the resonance curves A (w) induced by changes of 
the control parameters, leading to new nonlinear phenomena, have been stud- 
ied within the theory of algebraic curves - they occur in the neighbourhoods of 
singular points of A (w) [TT1 \12\ fT3 ] . 

In the present paper we study the exact fourth-order equation for internal 
motion. It turns out that the KBM method can be applied to yield equation 
defining the amplitude profile A (w; a, 6, . . .) implicitly, where a, 6, . . . are some 
parameters. This resonance curve is more complicated then in the case of ef- 
fective equation and hence more complicated metamorphoses are possible. The 
aim of the present paper is to explore these possibilities. 

The paper is organized as follows. In the next Section the exact 4th-order 
equation for the internal motion in non-dimensional form is presented. In Sec- 
tion 3 equation for the resonance curves A (w) is derived from the exact fourth- 
order equation for internal motion via the Krylov-Bogoliubov-Mitropolsky ap- 
proach. In Section 4 the theory of algebraic curves is used to compute singular 
points on the exact equation amplitude profiles - metamorphoses of amplitude 
profiles occur in neighbourhoods of such points. In Section 5 examples of an- 
alytical and numerical computations are presented for the 4th-order equation 
equation. Our results are summarized in the last Section. 

2 Exact equation for internal motion 

In new variables, x = Xi, y = X2 — X\, equations (TT]), @ can be written as: 

nix + vx + ax + V e (y) + R e (y) = f cos (ujt) 
m e (x + y)- V e (y) - R e (y) = 

where m = mi, m e = 7712, v = v\, a = ai, V e = V2, R e = R%- 

Adding equations ([3J we obtain important relation between variables x and 

y- 

Mx + v\x + a.\x + m e y = / cos (cot) , (4) 

where M = m + m e . 
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We can eliminate variable x in ^ to obtain the following exact equation for 
relative motion: 

[M^ + vf t + a) (py - V e (y) - R e (y)) + em e (vf t + a) y = Fcos (art) , (5) 

where F = m e w 2 /, /i = mm e /M and e = m e /M is a nondimensional parameter 
[10], see also Ref. [14] where separation of variables for a more general 
system of coupled equations was described. Equations ([5]) , (TJ| are equivalent to 
the initial equations (fTJ), (0. 

For small e we can reject the term proportional to e to obtain the approx- 
imate (effective) equation which can be integrated partly to yield the effective 
equation: 



mi - Ve (y) - Re (y) = 



■■ COS (uit + 8) 



V M2 (- 2 -#) 2 +" 2 

In what follows we shall assume 

R e {y) = ~a e y - -f e y 3 , V e (y) = -v e y + X e y 3 . 



(6) 



(7) 



This model was also investigated in [8] where limiting phase trajectories ap- 
proach was used. 

In this work we shall investigate the exact equation (0. We write Eqns. ([5]), 
in nondimensional form. Introducing nondimensional time r and rescaling 
variable y: 



T = tu, z = y s — , 



where: 



we get: 

where C is a linear operator: 

and nondimensional constants are given by: 

h = ^,b = h = — , n = - g = 

/JUI 7e Mlu to 



(It 2 



a, 



1 



(8) 
(9) 

-^Gil 2 cos (Or) (10) 
(11) 



/> K — 7 a 

a e V a. e to a e M 



(12) 
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3 Nonlinear resonances via Krylov-Bogoliubov- 
Mitropolsky method 

We apply the Krylov-Bogoliubov-Mitropolsky (KBM) perturbation approach 
[15] to the exact nonlinear fourth-order equation (|10p describing internal motion 
of the small mass. The equation (| 10[) is written in the following form: 

£(^f + ft 2 z) +e(a£z + g(z,zj) = 0, (13) 
where C is defined in (1111) and eg (z, z) is given by: 

eg = t(hf T -b(^f + z + z 3 ) - Q 2 Cz + n (H± +a)£§- ^Gtf cos (ilr) , 

(14) 

and 

6 2 - n 2 = ea. (15) 

Equation (fT3|) was prepared in such way that for e = the general solution, 

z{t) = Acos(Qr + <p) +Cexp (-§ (h - Va) t) + L>cxp f-± (h + VA^J r), 

A = H 2 —4a, with constant and arbitrary A, y>, C, D and, moreover, the solution 
for H . a > does not contain secular terms and z (t) — > A cos (fir + ip) for 

T > OO. 

We shall now look for 1 : 1 resonance using the KBM method. For small 
nonzero e the solution of Eqns. (fl"3"| ) - (fTB"j) and (J7J is sought in form: 



z = Acos(flT + ip) + ezi (A, tp, r) + . . . (16) 
with slowly varying amplitude and phase: 

— = eM 1 {A,tp) + ..., (17) 

^ = eN 1 (A,<p) + .... (18) 

Computing now derivatives of z from Eqns. ([T6|) . (fl7|) . (|18|l and substituting 
to Eqns. fTB")) - (US]), © and eliminating secular terms and demanding Mi = 0, 
iVi = we obtain the following equations for the amplitude and phase of steady 
states: 

—QA (pfl 2 -q) + jflA 3 (H - bail 2 + bfl 4 ) + ^Gfl 2 sin tp = 0, (19a) 
—A (n 4 - rn 2 + a) - f A 3 (a-n 2 + Hbn 4 ) + ^GQ 2 cos tp = 0. (19b) 

where p = h + H (k + 1), q = ah + H, r = hH + a (« + 1) + 1. 

Solving the system of equations (|19j) we get the implicit expressions for the 
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amplitude A (ft) and the phase (p (ft): 

A (ft) = ?g- (20a) 

k + 1 vc 2 + z? 2 

tan^(ft) = — , (20b) 

C* = OA (ft) (pft 2 - q) - fftA 3 (ft) (if - baft 2 + 6ft 4 ) (20c) 
D = A (CI) (ft 4 - rfl 2 + a) + §A 3 (ft) (a - ft 2 + Hbfl 4 ) (20d) 

Equation for the correcting term zi is of form: 

t + ft 2 ^i) = \HftA 3 sin ($) + \A 3 ((3 + a) ft 2 - 36) cos ($) , (21) 

where $ = 3ftr + 3<^(ft). Solving Eqn. (|2"Tj) and substituting to (fTB"|) we get 
finally: 

z = A (ft) cos (ftr + tp)- ^ A 3 (ft) bfl sin ($) + ^ A 3 (ft) cos ($) (22) 
where A (ft) , ip (ft) are given by Eqns.(|20p. 

4 General properties of the function A (Q) 

After introducing new variables, ft 2 = A, A 2 = Y, the equation (|20ap defining 
the amplitude profile reads 

L (X, Y; a, b, h, H, k, J) = XY (pX -q - §Y (H - abX + bX 2 )f , . 
+F (X 2 - rX + a + \Y (a - X + bHX 2 )f - JX 2 = 

where, as before, p = h + H (k + 1), g = ah + H, r = hH + a (k + 1) + 1. A new 

parameter J is a renormalized G, J = (j^ij G 2 . To obtain the corresponding 

expression for the effective equation ([6]) one can put G = so that J = j 2 

and then k = 0. 

Singular points of L (X, Y) are computed from equations [16] : 

L = 0, (24a) 
f§ = 0, (24b) 
# = 0. (24c) 

We can eliminate J from Eqns. (|24a[) . (|24b|) computing L - \X^ = 35 YA 
where 



A' = -186 2 # 2 A 4 r 2 - 24a 2 (n + 1)17 + 24bhH 2 X 3 Y - A8bHX 4 Y 
-18aAF 2 - 48aAy + 32HhX 3 + 32a 2 - 27b 2 X 5 Y 2 + 16H 2 X 
-32aA - 32a 2 (k + 1) X + 32a (k + 1) A 3 + 9H 2 XY 2 + 2AH 2 XY 
+48a 2 Y + 18a 2 F 2 + 48hbX 4 Y - 9a 2 b 2 X 3 Y 2 + 36ab 2 X 4 Y 2 
+16a 2 h 2 X - 16 (k + l) 2 H 2 X 3 + A8bH (k + 1) A 4 Y~ 
-32/iff (k + 1) A 3 - A8abhX 3 Y + 32A 3 - l6h 2 X 3 - 32A 4 + 24A 3 F 



(25) 
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to obtain simplified equations: 



K = 0, (26) 
f# = 0, (27) 

from which X, Y can be computed as functions of parameters a, 6, h, H, k and, 
finally, J can be computed from the last equation 

m = o. (28) 

Equations (|21)1) , ([2~T1) , (|2"8|) are still very complicated making analytical inves- 
tigation virtually impossible. We shall thus solve these equations numerically. 



5 Computational results 

In the present Section singular points of amplitude profiles - solutions of Eqns. 
(|26|) , ([27)) , (|28|) - are studied. More exactly, resonance curves with one singular 
point, two singular points on one curve, and with degenerate singular point are 
presented and metamorphoses of bifurcation diagrams are shown. 

5.1 Amplitude profiles with one singular point 

We have computed singular points for the following values of control parameters: 
k = 0.05, b = —0.001, H = 0.4, a = 5, h = 0.5 obtaining four physical solutions 
(i.e. with X > 0, Y > 0, J > 0): 



Table 1. 


X 


Y 


J 


n 


2.170 051 157 


1.357 255 661 


1.656 917 694 08 


1 


4.835 083103 


4.192 055 014 


1.036 434 99178 


2 


2.798 801078 


1.237140 868 


1.814 387 388 23 


3 


4.153 001386 


4. 680111331 


0.963 352 653 58 
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The first two solutions correspond to self-intersections, see Fig. [l] while the 
second pair represents isolated points. 

Metamorphoses of bifurcation diagrams which occur in the neighbourhood 
of self-intersections for the exact fourth-order equation are, for small k, quali- 
tatively similar to those studied for the case of 1 : 1 resonance in the effective 
equation in [TTJ [13] and are not shown here. 
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Figure 1: Amplitude profiles with singular points, k = 0.05, b = —0.001, H = 
0.4, a = 5, h = 0.5, J = 1.656 917 694 08 (left self-intersection, red curve, n = 1 
in Table 1), J = 1.036 434 99178 (right self-intersection, red curve, n = 2 in 
Table 1) and neighbouring curves (blue and green lines). 



5.2 Amplitude profiles with two singular points 

It is possible, tuning the parameters properly, to obtain amplitude profile with 
two singular points. 

Let, as before, k = 0.05, b = —0.001, a = 5, h = 0.5, H being arbitrary. 
We can compute, for some H , from (|2"6")l , (f2~T| X (H) 1 Y (H), then from Eqn. 
(|2"8|) we get Ji (H) and J2 (H) corresponding to two curves with one intersection 
each. The condition for a curve with two intersections is J\ — J2 for some H. 

To find this value of H we compute J%, J2 for two values of H, Hq = 
0.40, Hi = 0.55, and use linear extrapolation to compute H = H cr such that 
Ji (H cr ) = Ji (H cr ). In one step of this procedure we compute new value of 



ff(i+2) f rom known and H^ +1 \ J x (i+1) , J^ +1) solving linear 

system of equations for a^' l+1 \ 



where i = 0, 1, 2, Then the next value of H^ l+2 ^ is computed as 7j( l + 2 ) = 




^(i+l) a (i,i+l) + ^(i.-i+l) 



(29) 




. The convergence is quite fast, see Tables 2,3. 
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Table 2 


if) 








i 


2.170 051 157 


1.357 255 661 


1.656 917 694 


0.40 





2.222181 140 


1.452 883 358 


1.736 397 285 


0.55 


1 


2.250 807092 


1.503 767 806 


1.775 975 009 


0.611498 955 


2 


2.247348 900 


1.497 675 627 


1.771343 329 


0.604 644 433 


3 


2.247 330153 


1.497 642 563 


1.771318111 


0.604 606 880 


4 


2.247 330182 


1.497 642 613 


1.771318150 


0.604 606 937 
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Tabic 3 






tW 
J 2 


H d) 


i 


4.835 083103 


4. 192 055 014 


1.036 434 992 


0.40 





4.793 018 228 


2.753 798 647 


1.555 975 411 


0.55 


1 


4.620 272 267 


2.513 582 261 


1.798 606 879 


0.611498 955 


2 


4.641565 672 


2.538 879 342 


1.771466 643 


0.604 644 433 


3 


4.641680 711 


2.539 018 448 


1.771317923 


0.604 606 880 


4 


4.641680 536 


2.539 018 236 


1.771318150 


0.604 606 937 
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Figures below show convergence of two curves with one singular point to one 
curve with two singular points. 




0.0 0.2 0.4 06 0.8 1.0 12 1.4 1.6 1.3 2.0 22 2.4 2.6 2.8 3.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 22 2.4 26 2.8 3.0 

Figure 2: Convergence of amplitude profiles to critical amplitude profile with 
two singular points (red curve): convergence of curves from Table 2 (left figure, 
green curves), and from Table 3 (right figure, blue). 

Bifurcation diagrams computed for parameters in the neighbourhood of such 
resonance curve display presence of two singular points, see Figs. [3] 

In Figures[3]the parameters are k = 0.05, b = —0.001, a = 5, h = 0.5 in both 
cases and J = 1.745 481261, H = 0.6054 for Fig. 3a and J = 1.740 481261, 
H = 0.6034 for Fig. 36. 
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Figure 3: Bifurcation diagrams. The left figure corresponds the amplitude pro- 
file with singular point (left figure) with two cusps and to the nonsingular curve 
(right figure) with two gaps. 

These diagrams correspond to amplitude profiles shown in Fig. 2) 




H ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 

1.0 U 1.4 1.6 U 10 XL 2.4 2.6 U 3.0 X 



Figure 4: The amplitude profile with two singular points (left figure, red) and 
nonsingular curve (right figure, magenta). 



5.3 Merging two singular points into a single degenerate 
point 

It is possible, by smooth change of the parameters, to merge two singular points 
lying on the red curve in Figs. 2, 3. The resulting singular point is degenerate, 



9 




Figure 6: Metamorphosis of the bifurcation diagrams near amplitude profile 
with degenerate singular point. 

Solving Eqns. {3D]) for k = 0.05, J = 1.771318150 (these two param- 
eters correspond to the parameters of the critical red curve with two singu- 
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lar points) we get X = 3.113 090 974, Y = 2.087 620 813, h = 0.548 982 679, 
a = 4.538 990 962, 6 = -1.718 542 532 x 10~ 2 , H = 0.644 095 068, see Fig. [5j 

Bifurcation diagrams computed in the neighbourhood of the degenerate sin- 
gular point depend sensitively on small changes of parameters, Fig. [SJ 

6 Summary and discussion 

In this work we have studied dynamics of two coupled periodically driven os- 
cillators. The inner motion of this system has been described by the exact 
fourth-order equation ([5|) (or (fT0|) in nondimensional form). Applying the KBM 
method we have computed approximate resonance curves (amplitude profiles) 
A Although the KBM method is basically used for the second-order equa- 
tions we managed to apply it to the fourth-order equation since it was possible 
to eliminate secular terms and impose steady-state conditions. Dependence of 
the amplitude A on the forcing frequency Q is complex since A (51) is defined 
implicitly as an algebraic curve, L(X,Y) = 0, see Eqn. (f2l?|) . with polynomial 
function L depending on variables X = f2 , Y = A 2 and control parameters a, 
b, h, H, k, J in a complicated manner. 

In our previous paper we stressed that near singular points of algebraic 
curves, defining amplitude profiles, metamorphoses of bifurcation diagrams (and 
hence of dynamics) take place. In the present paper we have studied three cases 
of singular points of the resonance curves defined by Ean. ([2U)) : i) the case of one 
singular point fSection !5.ip . ii) the case of two singular points on one resonance 
curve (Section (15.210 . iii) the case of degenerate singular point (Section (|5.3t ). 
Indeed, dynamics of the system (ITU1) changes significantly in the neighbourhood 
of singular points of resonance curve L (X, Y) = 0. Singular points described in 
Section[5]are just the tip of the iceberg and thus we are going to study multitude 
of singular points of amplitude profiles ([23]) in our future work. 
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